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A Comparison of Radio Axis with Host Galaxy Plane Axis in Seyfert 

Galaxies 



Henrique R. Schmitt li2 ' 4 , Anne L. Kinney 1 , Thaisa Storchi-Bergmann 2 and Robert Antonucci 3 



ABSTRACT 

We use the radio axis as an indicator of the orientation of the obscuring torus in 
Seyfert galaxies, and analyze the difference between the position angles of extended 
radio structures and host galaxy major axis of Seyfert 1 and Seyfert 2 galaxies. We find 
that Seyfert l's are less likely to have extended radio structures along the host galaxy 
major axis, while Seyfert 2's have these structures distributed in most directions. We 
also find a zone of avoidance in the distribution of position angles; both Seyfert l's and 
Seyfert 2's seem to avoid close alignment between the radio axis and the host galaxy 
plane axis. These results are analyzed from the point of view of a model in which 
Seyfert l's have their obscuring torus axis aligned preferentially along the host galaxy 
disk axis, and Seyfert 2's have their torus axis laying at an intermediate angle between 
the galaxy disk and its axis. 
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1. Introduction 

Since the discovery by Antonucci & Miller 
(1985) of polarized broad emission lines in the 
nuclear spectrum of NGC1068, the interest in 
the Unified Model for Seyfert galaxies has grown 
considerably (see Antonucci 1993 and Urry & 
Padovani 1995 for a review of the model). This 
model assumes that Seyfert l's and Seyfert 2's 
both contain an active nucleus surrounded by a 
dusty torus, and that the angle through which 
the central engine is observed determines the 
classification of the object. Several pieces of ev- 
idence corroborate this model, like the obser- 
vation of polarized broad emission lines in sev- 
eral Seyfert 2 galaxies (Miller & Goodrich 1990), 
deficit of ionizing photons in Seyfert 2's (Wilson, 
Ward & Haniff 1988; Kinney et al. 1991) and the 
collimation of the nuclear radiation, observed as 
extended linear radio sources (Ulvestad & Wilson 
1989) and conically shaped Narrow Line Regions 
(Pogge 1989). 

In a recent paper Schmitt & Kinney (1996) 
studied the NLR shape of Seyfert galaxies, as 
predicted by the Unified Model, using archival 
high resolution HST [OIII] images. Their results 
show that the sizes of Seyfert l's NLR's are much 
smaller than the NLR's of Seyfert 2's would be 
if they were observed pole-on. The sample has 
not been selected by an isotropic property, but 
the similarity in radio and A5007A luminosities 
of Seyfert l's and Seyfert 2's suggests the com- 
parison may be rather fair anyway. This result 
apparently contradicts the Unified Model, from 
which we would expect both kinds of objects to 
have similar intrinsic NLR sizes. In order to solve 
this problem, Schmitt & Kinney (1996) propose 
that the Seyfert l's torus axis may be aligned 
preferentially along the host galaxy plane axis, 
while the Seyfert 2's torus axis may lie at an in- 
termediate angle between the galaxy plane axis 
and the galaxy plane. In this picture, because 
the amount of extended gas intercepted by the 
nuclear radiation is smaller perpendicular to the 



plane than at directions closer to it, the NLR 
appears more extended in Seyfert 2's than in 
Seyfert l's. 

The reason for the Seyfert 1 orientation distri- 
bution may be that even those objects which are 
seen from the polar nuclear torus direction, but 
nearly edge-on to the host galaxy are classified 
as Seyfert 2 because of obscuration by dust in 
the host galaxy plane. There is ample evidence 
that the broad line region in edge-on galaxies is 
reddened (de Zotti & Gaskell 1985) or entirely ex- 
tinguished (Keel 1980; Lawrence & Elvis 1982). 

In order to check the orientation of the torus 
axis in Seyfert l's and Seyfert 2's, in this pa- 
per we compare the position angle (PA) of the 
extended nuclear radio structures of Seyfert l's 
and Seyfert 2's with their host galaxy's major 
axis PA. We show that there is a lack of Seyfert 
l's with radio structures aligned along the host 
galaxy major axis (which would correspond to 
the case where the obscuring torus is nearly per- 
pendicular to the galaxy plane), while Seyfert 2's 
have radio structures distributed along almost all 
directions, consistent with the model proposed 
by Schmitt k Kinney (1996). 

2. The Sample and Measurements 

We selected from the literature a sample of 
46 Seyfert galaxies (15 Seyfert l's and 31 Seyfert 
2's) having high resolution radio maps and show- 
ing linear or slightly resolved radio structures, as 
defined by Ulvestad & Wilson (1984a). In Table 
1 we give the galaxy names, together with their 
activity type (Seyfert 1 or Seyfert 2), Morpho- 
logical Type, the mean numerical index (T) of 
stage along the Hubble sequence (as defined in de 
Vaucouleurs et al. 1991), distance, 6cm flux, log- 
arithm of the 6 cm power, PA of the radio struc- 
ture (PArad), PA of major axis (PAma), the dif- 
ference between PA^D and PAma (APA) and 
the inclination of the host galaxy. The distances 
were calculated using the galaxy's radial veloci- 
ties relative to the local group (de Vaucoulers et 
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al. 1991) and H =75 km s" 1 Mpc" 1 . The in- 
clinations were calculated from the axial ratios, 
assuming that the face-on galaxy is basically cir- 
cular in shape. 

The PA's of the extended radio emission were 
obtained from the references in Table 1, using 
their published values or measuring it on their ra- 
dio maps when the PA was not given explicitly. 
The host galaxies major axis PA's were mostly 
obtained from de Vaucoulers et al. (1991), with 
the remainder from references cited in Table 1. 
The PA's were checked by looking at the galax- 
ies on the Digitized Sky Survey Plates. For 
the galaxies without values for the major axis 
PA available in the literature and for the cases 
where the published value was wrong (NGC5929, 
MRK573, MCG-8-11-11), the PA was measured 
by fitting ellipses over the outer isophotes of the 
digitized galaxy image. 

3. Results 

Figure 1 shows the histogram of the differ- 
ences between the PA's of the radio axis and of 
the host galaxy major axis. We can see that 
there is a lack of small values for Seyfert l's, 
indicating that they are less likely to have ex- 
tended radio structures along the host galaxy 
major axis. In other words, the cases in which the 
radio axis lies in the plane of the host galaxy are 
rare in Seyfert l's, as expected. Meanwhile, the 
Seyfert 2's have APA values evenly distributed 
from 0° to 70°. However, it is interesting to 
note that both distributions have a lack of ob- 
jects with radio structures well aligned with the 
host galaxy minor axis. Applying a KS test to 
the data, we find that the hypothesis that the 
two groups of galaxies are drawn from the same 
parent population is rejected at the 99.0% level, 
or 99.1% when we exclude double nuclei galaxies 
(NGC1144, MRK110, MRK266 and MRK463E). 

We must ask if the Seyfert l's and 2's in the 
sample have similar intrinsic properties, or if the 
above results might be traceable to selection ef- 



fects. First we compare the logarithm of the 6cm 
radio power of the two groups, to see if they are 
similar in the two groups of galaxies, since the 
more luminous objects might be expected to have 
larger radio structures, which are easier to detect 
and measure. In Figure 2 we show the histogram 
of the logarithm of the 6cm radio power, where 
we can see that both groups have very similar 
distributions of radio powers, except for the high 
luminosity tail of Seyfert 2's. This high lumi- 
nosity tail had already been observed by Meurs 
k Wilson (1984) and Wilson k Ulvestad (1989). 
However, these objects are undistinguished in the 
position angle histograms, so they do not create 
the claimed effects. Note that the former refer- 
ence includes disk emission which may be signif- 
icant in low luminosity objects. 

Wilson k Tsvetanov (1994) have recently pro- 
posed that the obscuring torus axis may be 
aligned with the galaxy plane axis in late type 
galaxies, while the obscuring torus could have 
any orientation in early type systems. In order 
to check if our sample is biased towards Seyfert 
l's in late type galaxies and Seyfert 2's in early 
type galaxies, we show in Figure 3 a histogram 
of their morphological types, where the parame- 
ter T= —4 corresponds to Ellipticals, the earli- 
est morphological type, and T=5 corresponds to 
Sc's, the latest morphological type. The group 
distributions are almost equal, with the only dif- 
ference being the existence of some Seyfert 2's 
with T< —2 (as early as SO) without any cor- 
responding Seyfert l's in this region of the dia- 
gram. These objects cannot be distinguished in 
the position angle histograms, assuring us that 
the groups are not biased by different morpho- 
logical types. 

We also check for a trend for the obscuring 
torus axis to align along the minor axis in late 
type galaxies. In Figure 4 we show a plot of the 
morphological type versus APA. There is no sys- 
tematic trend for late type objects to have large 
APA values, which means that their radio struc- 
tures are not preferentially aligned with the host 
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galaxy minor axis. This result suggests that the 
trend found by Wilson & Tsvetanov (1994) was 
most likely due to the small number of objects in 
their sample. 

4. Discussion 

The results presented in Figure 1 can be in- 
terpreted from the point of view of the scheme 
proposed by Schmitt & Kinney (1996), where the 
Seyfert l's have their torus axis aligned preferen- 
tially close to the galaxy plane axis, while the 
Seyfert 2's have their torus axis aligned at an in- 
termediate angle between the galaxy plane and 
its axis. Here we describe a simple model, de- 
veloped in order to study the orientation of the 
projected radio structure (assumed to be aligned 
with the torus axis) due to its inclination rela- 
tive to the line of sight. In our model we as- 
sume a uniform distribution of orientation angles 
between the obscuring torus axis and the host 
galaxy plane axis. Then we predict what the 
observed values for APA would be, given that 
uniform distribution. 

Figure 5 shows a cartoon which represents our 
model. In this Figure, the galaxy disk lies on the 
X-Y plane and the elongated radio structure, rep- 
resented by the unit vector kj, makes an angle 
(3 with the Z-axis (galaxy plane axis). The elon- 
gated radio structure, is equally likely to be in 
any position along the circle "c", which is de- 
scribed by the azimuthal angle 9, measured from 
the X to the Y-axis. The galaxy is inclined rel- 
ative to the line of sight (Z'-axis) by an angle 
"i", such that the inclined galaxy major axis is 
coincident with the X-axis. 

In this model the angle (3, between the radio 
axis and the galaxy plane axis assumes values in 
the range 1° < (3 <90° and is varied in steps of 
1°. For each value of (3 we generate one thousand 
vectors equally spaced along one half the circle 
"c" . Due to the symmetry of the problem, the 
angle 9 is varied only in the range -90° < 6 <90°. 
Also, due to the fact that the circle "c" describes 



different perimeters for different values of (3, the 
results were weighted for each value of (3, by 
the area of a ring of 1° described by the circle 
"c" on the surface of a sphere, relative to the 
whole sphere area. This is necessary in order to 
avoid an oversampling in the number of vectors 
for smaller [3 values, relative to larger ones. The 
galaxy inclination "i" is then varied in the range 
15° <i<60° in steps of 1°. We chose this range 
because for i<15° it is difficult to measure the 
galaxy inclination angle and it is considered to 
be face on, while for i>60°, there is only a small 
number of galaxies with such inclination in our 
sample. The small number of objects with i> 60° 
suggests an inclination dependent selection bias. 
We have also tested the model for smaller and 
larger values of "i" , but the final result does not 
change. 

The projected components of each individual 
radio jet vector (for every 9, (3 and i value), in 
the directions X, Y' and Z', can be calculated by 
the following relations: 

Px = cos 9 sin (3 
Pyi = cos i sin 9 sin (3 + sin i cos (3 

Pz> = cos i cos (3 — sin i sin 9 sin (3 

The observed difference between the elongated 
radio structure PA and the major axis PA (APA) 
is given by: 

APA = arctan(Py//P x ) 

We divide APA in 10° bins and count the num- 
ber of vectors with projected APA inside each 
bin, for a given (3 value. This number, normal- 
ized to the total number of vectors, can be con- 
sidered as the probability of an elongated radio 
structure, which makes an angle (3 with the host 
galaxy axis, being seen at a projected angle APA 
from the galaxy major axis, considering inclina- 
tion effects. 
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The model also includes some constraints. 
First, a galaxy is only considered to be a Seyfert 
1 if the angle between Z' and k, is less than 30°, 
corresponding to observing into the opening of 
the torus; otherwise it is a Seyfert 2. This num- 
ber was obtained from Osterbrock & Shaw (1988) 
and corresponds to half the torus opening angle. 
Second, if the projected length (r) of the unit vec- 
tor kj in the X-Y' plane is less than a given value, 
it is assumed to be unresolved and is not consid- 
ered when we count the number of vectors that 
fall inside different APA bins. This corresponds 
to the case where the elongated radio structure 
is observed pole on, and is therefore not observed 
as elongated. Since we do not know the intrinsic 
length of the elongated radio structure, we note 
that they usually have sizes of the order of l"-2" 
and that the VLA resolution at 6cm in A config- 
uration is «0.4". We thus study two cases: r>0.2 
and r>0.4. Notice also, that Seyfert 2's are not 
effected by the projected size restriction, because 
the vectors with angle between Z' and kj larger 
than 30° will have r>0.5. 

Figure 6 shows the model results. The in- 
dividual lines are cumulative histograms, repre- 
senting the sum over all (3 angles smaller than 
or equal to the value indicated at the right and 
below each line. The vertical axis of these plots 
were normalized relative to the APA bin with the 
largest number of vectors on the f3 =90° bin. The 
histogram can be considered as the cumulative 
probability of an elongated radio structure, that 
makes an angle smaller than or equal to (3, being 
seen at a projected angle APA from the galaxy 
major axis. In Figure 6a we show the results for 
Seyfert l's with projected length (r) larger than 
0.4. Figure 6b shows the results for Seyfert l's 
with projected length larger than 0.2 and figure 
6c shows the results for Seyfert 2's. The APA's 
distribution for Seyfert l's are very similar, inde- 
pendent of the "r" value. For small (3 values there 
is a larger probability of observing an elongated 
radio structure closer to the host galaxy minor 
axis (large APA). For larger (3 values, the proba- 



bility of observing the elongated radio structure 
closer to the major axis increases. When we con- 
sider the distribution for all (3 angles together 
(the 90° line) , the probability of finding a vector 
at any APA bin is approximately the same. The 
distribution of APA's for Seyfert 2's is similar 
to that of Seyfert l's. However, the probability 
of finding small APA's, which corresponds to an 
elongated radio structure close to the host galaxy 
major axis, only increases when we go to larger 
(3 values relative to those of Seyfert l's. 

Comparing the model results from Figure 6 
with the observed APA's in Figure 1, we see that 
the lack of Seyfert 1 galaxies with small APA val- 
ues can only be understood if there are very few 
Seyfert 1 galaxies in which the angle between the 
elongated radio structure and the galaxy plane 
axis (0 in the model) is large {(3 > 30° - 40°). 
Similarly, when we compare the results for the 
Seyfert 2 models with the observed values, the 
distribution of APA values can be understood if 
the Seyfert 2 galaxies can take on any value for 
the angle between the elongated radio structure 
and the galaxy plane axis. 

These results support the scheme proposed by 
Schmitt & Kinney (1996). However, that scheme 
cannot explain the lack of both Seyfert l's and 
Seyfert 2's with large APA values, which would 
correspond to elongated radio structures aligned 
with the host galaxy minor axis. This lack of 
large APA values suggests a physical zone of 
avoidance, already observed by Ulvestad & Wil- 
son (1984b) on a smaller sample, where for some 
reason the radio axis is not ever closely aligned 
intrinsically with the galaxy plane axis. If we 
consider a flat distribution of APA's, we calcu- 
late that for our sample of 46 galaxies, we would 
expect 10 galaxies with APA> 70°. Considering 
a sample of only 31 galaxies, which corresponds 
to the number of Seyfert 2's in our sample, we 
would expect 7 galaxies with APA> 70°. From 
Poisson statistics we calculate that the a posteri- 
ori probability of observing only one galaxy with 
APA> 70°, when the expected number is 10, is 
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10 4 - 78 , or 10 3 47 for the case when 7 galaxies 
are expected. 

There may be some effect, which we are not 
taking into account, that makes large APA un- 
observable. One possible solution for the lack 
of Seyfert l's with elongated radio structures 
aligned with the host galaxy minor axis could be 
that when the obscuring torus is closely aligned 
with the host galaxy axis, an elongated radio 
structure is not produced. Due to the low HI den- 
sity in the central region of spiral galaxies (Rots 
1975, Begeman 1988), the radio "jet" would not 
interact with an enough large quantity of mat- 
ter and consequently would not radiate enough 
to be detected. Alternatively, these nuclear disks 
in the same plane as the host galaxy may have 
fewer dissipation/fueling mechanisms. Yet an- 
other possibility is that the nuclear axes reflect a 
past triggering merger, and that the approaches 
of companions are statistically anisotropic. To 
incorporate such effects, we show in Figure 7a 
and b (r>0.4 and r>0.2 respectively) the mod- 
els for Seyfert l's with f3 <10° excluded. In this 
case the probability of observing an elongated ra- 
dio structure with APA>70° is smaller than for 
values of APA< 70°. This makes the observed 
and modeled APA distributions look much more 
alike, although we would still expect to detect 
some objects with APA>70°. While this solu- 
tion can perhaps solve the problem for Seyfert 
l's, it does not suffice for Seyfert 2's, where the 
statistics are better. A simple visual inspection 
of Figure 6c shows that the exclusion of all vec- 
tors with (3 < 10° does not change the distribu- 
tion of APA's significantly. Even if we exclude all 
vectors with (3 <30°, the probability of observing 
a Seyfert 2 with APA>70° would be smaller, but 
it would not explain the deficit that we observe. 

The lack of Seyfert 1 galaxies with elongated 
radio structures aligned with the galaxy major 
axis is consistent with results obtained by Keel 
(1980) and Maiolino & Rieke (1995). These au- 
thors showed that there is a deficiency of edge-on 
Seyfert 1 galaxies and that Seyferts 1.8 and 1.9 



are more likely to be in edge-on galaxies. These 
results suggest that in the case of edge on objects 
we can be observing the nuclear region directly, 
through the polar region of the nuclear torus, but 
not seeing the broad line region due to shadowing 
by gas and dust in the galaxy disk. 

On the other hand, our results differ from 
those from Colbert et al. (1996) and Baum et al. 
(1993). These authors found that the large scale 
(^1 kpc) radio structure of Seyferts are prefer- 
entially aligned with the host galaxy minor axis. 
Baum et al. (1993) have also compared the small 
scale radio structures with the large scale radio 
structures and found that their PA's are differ- 
ent. Colbert et al. (1996) suggest that the small 
scale radio jets are possibly diverted by dense 
molecular clouds, like the scenario proposed for 
NGC1068 by Gallimore, Baum & O'Dea (1996), 
and then generate the large scale radio structures 
that we see. 

5. Summary and Conclusions 

We have shown in this Paper that there is a 
lack of Seyfert 1 galaxies with radio structures 
aligned with the host galaxy major axis, while 
for Seyfert 2's the radio structures are oriented 
along almost any direction in the galaxy. Both 
groups also show a deficiency of objects with 
elongated radio structures closely aligned with 
the host galaxy plane axis. 

We developed a model to calculate the distri- 
bution of APA, the difference between the ori- 
entation of the radio axis and host galaxy major 
axis, based on the assumption that the angle be- 
tween the radio axis and the host galaxy plane 
axis ((3) is uniformly distributed between 0° and 
90°, taking into account the effect of the galaxy 
inclination and resolution of the elongated radio 
structure. From the comparison of the observed 
APA distribution of Seyfert l's with the distri- 
bution predicted by the model, we can explain 
the small number of objects with small APA's 
only if their torus axis lies closer than ss 30° 
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to the host galaxy axis. The lack of Seyfert l's 
with APA> 70° can be partially explained if we 
assume that the elongated radio structures are 
not often produced closer than rs 10° to the host 
galaxy plane axis, and several possible explana- 
tions for this are suggested. For the Seyfert 2's, 
the observed distribution can be explained if the 
torus axis assumes any angle relative to the host 
galaxy plane, with the exception of the region 
with APA> 70°. 

These results, together with those from Keel 
(1980) and Maiolino k Rieke (1995), showing 
that Seyfert l's are less likely to be found in edge- 
on galaxies, as well as the results from Lawrence 
k Elvis (1982) and deZotti k Gaskell (1985), 
showing that the broad line region in edge-on 
galaxies is reddened, are in good agreement with 
the model proposed by Schmitt k Kinney (1996). 
The paucity of objects showing close alignment 
between radio axis and galaxy plane axis remains 
unexplained. 
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Fig. 1. — Histogram of the difference between 
the radio and the host galaxy major axis position 
angles of Seyfert l's (dashed line) and Seyfert 2's 
(solid line). 

Fig. 2. — The distribution of the logarithm of the 
6cm radio power. Lines as in Figure 1. 

Fig. 3. — The distribution of the galaxies mor- 
phological types (T). Lines as in Figure 1. T= 
—4 corresponds to morphological type E, T= —2 
corresponds to SO, T=0 corresponds to SO/a, 
T=2 corresponds to Sab and T=4 corresponds 
to Sbc. 

Fig. 4. — Comparison between the radio axis 
and major axis position angle differences and the 
host galaxy Morphological Type. Open circles 
are Seyfert l's and filled squares are Seyfert 2's. 

Fig. 5. — This cartoon shows our radio "jet" 
model. The host galaxy disk lies in the X-Y 
plane, the radio "jet" is represented by the vec- 
tor kj, which makes an angle (3 with the Z-axis 
(galaxy plane axis) and is likely to be in any di- 
rection along the circle "c" , which is described by 
the angle 6. The angle "i" describes the galaxy 
inclination relative to the line of sight. 

Fig. 6. — Results of the models. The lines are 
the cumulative sum of all the vectors with (3 less 
and equal the value indicated bellow the line, a) 
Syl's with r>0.4; b) Syl's with r>0.4; c) Sy2's 

Fig. 7. — Results of the Syl models for the case 
when we exclude the cases with (5 less and equal 
to 10°. a) r>0.4; b) r>0.2. 
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